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We propose multilayer cylindrical invisibility cloaks that are optimized for 
oblique incidences through combination of analytic formalism of scattering 
and genetic optimization. We show that by using only four layers of homo- 
geneous and anisotropic metamaterials without large values of constitutive 
parameters, the scattering for oblique incidences can be reduced by two 
orders. Although the optimization is done at a single incident angle, the cloak 
provides reduced scattering over a large range of incident angles. © 2010 
Optical Society of America 

A lot of efforts have been made recently on designing and implementing metamaterial 
invisibility cloaks using transformation methods [1-9]. A perfect invisibility cloak created 
from transforming empty electromagnetic space [1] has been proven theoretically to be perfect 
in cloaking an arbitrary object [5-7], while the only electromagnetic mechanism so far to 
detect a perfect cloak within its working band needs extreme conditions [10]. However, the 
practical development of metamaterial invisibility cloaking is still far from real application 
at this stage, since the rigorous requirements of an ideal cloak, such as continuously varying 
parameters which approach extreme values in some regions, are challenging for practical 
fabrication. Although a simplified cylindrical cloak that was implemented previously [3] 
can moderately reduce the scattering, it is still inherently visible [11], requiring further 
improvements on its performance. 

Moreover, previous studies on the cylindrical invisibility cloak, which is often preferred in 
most analysis and experiments because of its relatively simple geometry, were mostly limited 
to an in-plane or purely two-dimensional (2D) case where the wave must propagate in a 2D 
plane perpendicular to the cylindrical axis, which has severely hindered the practicability of 
cylindrical cloaking. On the other hand, recently, there have been theoretical suggestions to 
overcome the drawbacks of continuous and extreme parameters by optimizing the structure of 
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Fig. 1. Configuration of a 4-layer cloak when a vertically polarized plane wave is obliquely 
incident with the incident angle of a. 

a few layers of homogeneous and anisotropic metamaterials to construct a practical cloak [12, 
13]. However, similar to previous studies, only normal incidence is considered and therefore 
the results would only work for a purely 2D geometry. The extension of a 2D design to a 3D 
design is one step further towards the practically usable realization of invisibility cloaking. 

We present an invisibility cloak design with only 4 layers of metamaterials that can work 
at a wide range of incident angles. Our work can be treated as the extension and supplement 
of the design of purely 2D models into more physical and practical 3D models. We use 
genetic optimization to globally search an optimized set of parameters. The optimization of 
a multilayer cloak is done with oblique incidence of 30 degrees. Our method can be used as 
a guidance for future implementation of metamaterial invisibility cloaks. 

The geometry of the multilayer cylindrical cloak model is shown in Fig. 1, where a cylin- 
der of perfect electric conductor (PEC) with radius Rm+i is cloaked by the multilayered 
cylindrical shell with innermost radius Rm+i and outermost radius R\. In between R\ and 
Rm+i are M layers marked by m (1 < m < M) with different outer radius varying from Ri 
to Rm- The constitutive parameters (e p , e^, e 2 , /i^ and jj z ) of each layer are assumed to 
be homogeneous. The region of r > Ri is free space. The incident electromagnetic wave is 
incident on the multilayer cylindrical cloak with an angle of a with respect to the xy plane. 
Without loss of generality, we consider vertical polarization as shown in Fig. 1. Therefore, 
the unit incident wave takes the form of Ei = (— xsina + z cos a )e lhsma+lfacosa . 

For oblique incidence, the difficulty here is that there is no closed-form solution to the 
wave equation inside each layer. We have provided a method to calculate the scattering 
from a general cylindrical cloak under oblique incidence in [14] based on the state-variable 
approach [15]. We apply this method in each layer to calculate each corresponding state 
propagator matrix, respectively, and their product will be the final state propagator matrix 
of the entire multilayer structure. If we divide each layer into N sublayers (N 1), the state 
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propagator equation will be 



V(R M+ i) = 

nf= N (M-i)N+i(t + a p^(pj)) 

■v(Ri) (i) 

where the four-dimensional vector V = [E z H z H^] 7 * is the state vector. After obtaining 
the state propagator matrix, we can solve the field distribution over the entire space [14]. Now 
by applying the genetic optimization [12], a multilayer cloak applied at oblique incidence is 
obtained. 

For comparison with the performance of simplified cloaks at oblique incidence, a cylindrical 
structure with the same size as in [14] is considered, i.e. R± = 1.5Ao = 2.08/?,a/+i- We set 
the number of layers M to be 4 where each layer has the same thickness. Another condition 
is that we force all relative constitutive parameters in each layer to be less than 10. The 
incident angle of the incident wave is chosen to be a = 30°. In the genetic algorithm, the 
far-field scattering efficiency Q sca [7], or the scattering cross section normalized by the unit 
geometrical cross section of 2Rm+i, is chosen as the function to be minimized. 

We first plot the scattering from the bare PEC core (radius: 0.721Ao) as the reference, 
where the E z field in the xy plane is shown in Fig. 2(a). It is shown that the scattering is 
quite large. The scattering efficiency is 1.175 in this case. Note that although the incident 
wave is strictly vertically polarized, the reflected wave is generally composed of both vertical 
and horizontal polarizations, where only vertical polarization has E z component. The next 
task is to minimize Q sca by optimization. 

We have two methods of optimization. One is to control e z , fi p and ^ in each layer 
(the previous simplified cloak for normal incidence only need these three parameters), while 
keeping e p , and \i z as constant of 1. The other is to control all six constitutive parameters in 
each layer. The performances of these two optimized 4-layer cloaks are shown in Fig. 2(b-c). 
Their relative constitutive parameters are summarized in Table 1. It can be seen in Fig. 2(b) 
that by controlling only three parameters in each layer, the total scattering has reduced 
much when compared to the bare PEC core in Fig. 2(a). The scattering efficiency Q SC a is 
reduced from 1.175 to 0.498. After optimizing all six parameters in each layer in the cloak, 
the scattering efficiency Q sca is 0.013, close to "perfect invisibility" as shown in Fig. 2(c). It 
is worth mentioning that although the 4-layer cloak in Fig. 2(c) is achieved by optimization 
at the single incident angle of 30°, it is still valid for other incident angles. In Fig. 3, we can 
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Fig. 2. (Color online) field distribution in plane due to scattering from (a) a bare 
PEC cylinder, (b) a 4-layer cloak with three parameters (e z , fi p and /i^) optimized, and (c) 
a 4-layer cloak with all six parameters optimized, when a vertically polarized plane wave 
is obliquely incident from left to right with the incident angle of 30°. From (a) to (c), the 
scattering efficiency Q sca is 1.175, 0.498 and 0.013, respectively. 
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Fig. 3. (Color online) Dependence of scattering efficiency Q sca on the incident angle for the 
4-layer cloak achieved by optimizing all six parameters in each layer within the cloak. The 
dotted line corresponding to the scattering efficiency of the bare PEC core is shown for 
comparison. 

see that this cloak has satisfactory performance over a large range of incident angles. 

In principle, we can ascribe this low scattering phenomenon to the destructive interference 
between the waves from the first reflection on the outermost boundary and the waves from 
the total transmission from inside the cloak to outside [12]. However, the scenario of oblique 
incidence is more complicated than the case of normal incidence, because both polarizations 
need to be optimized at the same time. For normal incidence, since vertically polarized waves 
and horizontally polarized waves are decoupled, reflection and transmission at individual 
boundaries between adjacent layers will affect one polarization without influence from the 
other polarization. Therefore it is relatively easier to achieve destructive interference outside 
the cloak. When the incident angle is nonzero, vertically polarized and horizontally polarized 
waves become coupled. Thus there is a need for more degrees of freedom to optimize the 
structure to achieve destructive interference for both polarizations. This justifies a much 
better performance of the 4-layer cloak with all constitutive parameters optimized in each 
layer (Fig. 2(c)), as oppose to just 3 parameters (Fig. 2(b)). 

In conclusion, we have proposed a method to achieve multi-layered cylindrical cloaks work- 
ing at oblique incidence by combining the analytic algorithm of scattering and genetic opti- 
mization. We show that by using only four layers of homogeneous and anisotropic metama- 
terials with constitutive parameters of finite values, the scattering from the cloak at oblique 
incidence can be reduced by two orders. Although the optimization is performed only at a 
single incident angle, the performance of the cloak is satisfactory within a large range of 
incident angles. We anticipate that the combination of analytic algorithm of scattering and 



5 



I 


Optimizing three parameters 


layer 


Hp 


H<f> 


e z 


e p e <t> 


Hz 


1 


0.962 


1.423 


1.154 


1 1 


1 


2 


1.128 


0.160 


0.598 


1 1 


1 


3 


0.792 


2.754 


1.321 


1 1 


1 


4 


0.665 


1.025 


2.625 


1 1 


1 


II 


Optimizing six parameters 


layer 


H P 


H<t> 


e z 


e P 


Hz 


1 


0.550 


2.320 


1.725 


0.438 1.781 


3.409 


2 


0.283 


3.584 


1.657 


0.262 3.953 


0.037 


3 


0.224 


8.204 


0.614 


0.189 8.674 


1.754 


4 


0.057 


9.994 


0.015 


1.540 9.069 


0.030 



Table 1. The relative constitutive parameters for the optimized 4- layer cloak by (I) optimizing 
three parameters in each layer and (II) optimizing six parameters in each layer. 

genetic optimization as a powerful tool will be used in more applications in the practical 
implementation of invisibility cloaks in the future. 
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